arXiv: 1501.01175v3 [cond-mat.str-el] 4 Mar 2015 


Probing the Fermi Surface of 3 D Dirac Semimetal Cd3As2 through de Haas-van 

Alphen Technique 

A. Pariari, P. Dutta and P. Mandal 
Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta 700 064, India 

(Dated: March 5, 2015) 

We have observed Shubnikov-de Haas and de Haas-van Alphen effects in the single crystals of 
three dimensional Dirac semimetal Cd 3 As 2 upto to 50 K, traceable at field as low as 2 T and 1 
T, respectively. The values of Fermi wave vector, Fermi velocity, and effective cyclotron mass of 
charge carrier, calculated from both the techniques, are close to each other and match well with 
the earlier reports. However, the de Haas-van Alphen effect clearly reflects the existence of two 
different Fermi surface cross-sections along certain direction and a non-trivial Berry’s phase which 
is the signature of 3D Dirac Fermion in Cd 3 As 2 . 


Now a days, topology-dependent electronic properties 
of solids are the subject of considerable research interest 
in condensed matter physics. Dirac materials such as 
graphene [I| , topological insulator [3, 0| and the parent 
compounds of iron based superconductors whose 
excitations obey a relativistic Dirac-like equation have 
been widely studied in recent years. Three dimensional 
(3D) Dirac semimetal is one of the candidates of topol¬ 
ogy protected band structure in which semimetal bulk 
has linearly dispersing excitation and surface state is the 
topology protected Fermi arc [1-Q. In these systems, 
the pair of 3D Dirac points is protected by crystal 
symmetry [iSIl. From Dirac semimetal, it is possible to 
obtain a Weyl semimetal or topological insulating phase 
by breaking time reversal or inversion symmetry [l^ . 
One can also obtain topological superconducting phase 
by doping charge carrier in the system [l^ . Recent 
experiments have demonstrated the existence of a 3D 
Dirac semimetal phase in NaaBi and Cd 3 As 2 [IMl- 
Angle-resolved photoemission spectroscopy (ARPES) 
measurements on Cd 3 As 2 single crystals have shown 
that two almost identical ellipsoidal Fermi surfaces 
are located on the opposite sides of the Brillouin zone 
center point along the kz direction BUS- However, 
these two are expected to merge into a single ellipsoidal 
contour around F point beyond Lifshitz transition 
from STM [I^. The Lifshitz transition occurs when 
the two Fermi surfaces touche with each other. Band 
structure calculations suggest that this transition should 
occur very close to the Fermi energy~133 meV [ 11 - 
Analysis of Shubnikov-de Haas (SdH) oscillations also 
suggests the existence of one or equivalent Fermi surface 
cross-sections corresponds to the single frequency of 
oscillation Thus, probing the Fermi surface of 

Cd 3 As 2 by different techniques may provide a crucial 
cross reference on its novel quantum properties. 

Shubnikov-de Haas (SdH) oscillations, which have 
been observed [isl-l^ in Cd 3 As 2 , originate due to the 
oscillations in the scattering rate of charge carriers. 
On the other hand, the oscillations in free energy in 
presence of magnetic field cause de Haas van Alphen 
(dHvA) effect, i.e., the oscillations in magnetization. 


Though the oscillations in both the cases are driven 
by the same physical parameter, dHvA offers a more 
accurate picture of the Fermi surface as it is not 
sensitive to quantum interference effects or noise from 
the electrical contacts. However, the investigation of 
Fermi surface in this 3D Dirac semimetal through dHvA 
is yet to be done. In this Letter, along with SdH 
oscillation and high non-saturating linear magnetoresis¬ 
tance (MR), we report the observation of two different 
Fermi surface cross-sections from the dHvA effect 
which eventually establishes itself as a better probe to 
study the Fermi surface properties. Furthermore, this 
non-transport dHvA effect also confirms the nontrivial 
Berry’s phase (signature of 3D Dirac fermion) in Cd 3 As 2 . 

Single crystals of Cd 3 As 2 were synthesized by chemical 
vapor transport technique. First, we prepare polycrys¬ 
talline Cd 3 As 2 samples by heating the stoichiometric 
mixture of high purity Cd pieces and As powder at 
500°C for 8 h and at 850°C for 24 h in a vacuum sealed 
quartz tube of 15 cm long and 16 mm diameter. The 
product was then heated again at 550° C for another 
48 h for homogenization. Finally, the quartz tube was 
placed in a gradient furnace and heated for 48 h. During 
heating, the end of the quartz tube which contains the 
compound was maintained at 690°C while the other 
end was kept at 600°C. The furnace was then cooled 
slowly to room temperature. Several small size shiny 
plate-like crystals formed at the cold end of the tube 
were mechanically extracted for transport and magnetic 
measurements. Powder X-ray diffraction of crashed 
single crystals shows that these crystals have Ri/acd 
space group and contain no impurity pha ses. The details 
are given in supplementary material [23|. The electrical 
and thermal transport measurements were carried out in 
a 9 T physical property measurement system (Quantum 
Design). The magnetization was measured using a 
SQUID VSM (Quantum Design). 

Figure 1(a) shows the temperature dependence of 
the resistivity pxx of a Cd 3 As 2 single crystal, pxx is 
metallic over the whole range of temperature. pxx 
exhibits approximately ® dependence above 100 
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FIG. 1. (a) Temperature dependence of resistivity of Cd 3 As 2 
crystal, (b) Magnetic field dependence of magnetoresistance 
(MR) at various temperatures when field along [100] direc¬ 
tion.(Inset temperature dependence of MR at field 9 T.) 


K and becomes almost T independent for T<20 K 
which reveals a residual resistivity pxxiO) ^ 50 pfl 
cm. Both the value and the nature of temperature 
dependence of resistivity are similar to earlier reports 
(Ta [ 2 ^. Fig. 1(b) displays the field dependence of the 
magnetoresistance at different temperatures. MR is 
defined as [pxxiB) - Pxxi^)]/Pxx{^)- Here, current is 
along [012] direction and the applied magnetic field is 
parallel to [100] direction. Even at room temperature 
and 9 T magnetic field, MR shows no sign of saturation 
and its value is as high as 263%. Except at low field, 
MR is approximately linear. Such a large linear MR at 
room temperature is quite unusual and needs further 
theoretical study for understanding the physical origin. 
With decreasing temperature, MR increases rapidly and 
it reaches ~1600% at 2 K and 9 T. In the inset of Fig. 
1(b), we have shown the temperature dependence of 
MR at 9 T in log-log scale. It is clear from the figure 
that, unlike conventional metal, MR shows a saturation 



FIG. 2. (a) The oscillatory component AR^ja, of MR(i3 along 
[100] direction) [after subtracting a smooth background] as a 
function of 1/B at various temperature and the inset shows 
the oscillation frequency after fast Fourier transform, (b) 
ARa;a; for magnetic field along [021] and the inset shows the 
corresponding frequency. 


like behavior in the low temperature region below 30 K. 
Similar behavior has also been observed in other crystals. 

ARxx, obtained after subtracting a smooth back¬ 
ground from MR, is plotted in Figs. 2(a) and (b) as a 
function of 1/B for fields along [100] and [021], respec¬ 
tively. The presence of SdH oscillations can be traced 
down to a field of ~2 T and up to 50 K. With increasing 
temperature, the amplitude of oscillation decreases 
rapidly. The fast Fourier transform spectra of AR^x 
versus l/B curve shows a single oscillation frequency 
at around 56 T [Inset of Fig. 2(a)] and 53 T [Inset of 
Fig. 2(b)] for held along [100] and [021], respectively. 
According to the Onsager relation F={(j)o/2TT^)Ap, the 
cross-sectional areas of Fermi surface normal to the held 
are Ai;’=5.32x 10“^ and 5.04 x 10“^A respectively. 
Assuming a circular cross-section, a small Fermi momen¬ 
tum kp ~0.04 A is calculated. Both the values and 
presence of single frequency are in good agreement with 
recent ARPES [H, [l^ and STM jlSl ] results. Where 
in both the cases single frequency is expected for two 
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identical ellipsoids or single ellipsoidal contour, below 
and above the Lifshitz transition, respectively. 

The magnetization of Cd 3 As 2 crystal has been mea¬ 
sured at different temperatures with B parallel to [021] 
direction, which displays a very clear dHvA effect [Fig. 
3(a)]. We have been able to observe the oscillations 
down to 1 T and up to 50 K. This indicates that the 
quantization of electron orbit does not get blurred 
by collisions with phonon or any impurity due to the 
high mobility of the charge carriers. In Fig. 3(b), the 
magnetic susceptibility Ay {=dM/dB) versus \/B plot 
shows how the oscillation amplitude decreases with 
increasing temperature. Inset of Fig. 3(b) shows the 
oscillation frequency obtained by the Fourier transfor¬ 
mation of dM/dB curve. Unlike magnetoresistance 
data, we observe that there are two distinct oscillation 
frequencies: one at 46 T, and the other at 53 T. This 
reveals two cross-sectional areas of the Fermi surface 

perpendicular to B, 4.39 x 10“^ and 5.04 x 10“^ 
° —2 

A respectively. To understand the phenomenon, 
we have done the magnetization measurements with 
field along [100] direction [l^. But dHvA oscillation 
along this direction gives only one frequency which 
is close to the frequency found from SDH oscillation 
along that direction. Though the dHvA oscillation 
along [100] direction establishes the equivalence between 
two ellipsoidal Fermi surfaces found in ARPES 0 
or single ellipsoidal contour after Lifshitz transition, 
the magnetization measurement along [02l] direction 
indicates two Fermi surface cross-sections correspond to 
the two frequencies. Also, we have done magnetization 
measurements in [100] plane, by applying magnetic field 
making an angle 9 with [021] direction. It is found that 
the peak amplitude at 53 T diminishes with increasing 
9 and disappears at 0 = 90° i.e., along [012] direction, 
giving single frequency of oscillation at 46.5 T (^ . 
The above study clearly shows that two equivalent 
cross-sections from two ellipsoidal Fermi surfaces super¬ 
pose with each other along certain direction leading to 
two frequencies of oscillation in dHvA effect. Another 
possible explanation for the occurrence of these two 
frequencies may be the contribution from surface and 
bulk states as predicted theoretically [IJ . But the single 
frequency along [100] and [012] directions and sample 
thickness ^0.4 mm which seems to be much larger 
than the critical thickness, eliminate the possibility of 
quantum oscillation from surface state. 

In Figs. 4(a) and (b), the integer index n, corre¬ 
sponding to the peak in SdH and n-|-l/4 corresponding 
to peak in dHvA oscillation, is assigned in such a way 
that a linear extrapolation of the straight line yields 
an intercept closest to zero index. Here we have taken 
into account the fact that Ay = dM/dB is advance 
in phase by tt/2 with respect to AM (oscillation in 
magnetization after subtracting a smooth polynomial 
back ground) [ 2 ^. In order to determine the quantum 




FIG. 3. (a) Showing magnetic field(along [021] direction) de¬ 
pendence of diamagnetic moment at temperature 2 K and 
20 K,(b) Oscillating part of susceptibility Ay= d{M)/dB vs. 
l/B plot (de Haas van Alphen (dHvA) effect) and inset shows 
oscillation frequency after fast Fourier transform found from 
dHvA effect. 


limit, i.e., the value of magnetic field required to achieve 
the lowest Landau level (n=0), we have plotted n vs l/B 
in Figs. 4(a) and (b). From the linear extrapolation of 
index plot of SdH and dHvA effect, it is found that the 
system reaches quantum limit at~ 50-60 T. The value 
of critical field for the quantum limit is very close to 
that observed recently [^, [^ . When a closed cyclotron 
orbit is quantized under an external magnetic field B, 
then according to the Lifshitz-Onsager quantization rule 
Api^/eB) = 2TT{n + \/2 — (3 +d) Here, 27r/3 is the Berry’s 
phase where /3 can take values from 0 to 1/2 (0 for 
parabolic dispersion as in the case of conventional metals 
and 1/2 for the Dirac system), and 2ti5 is an additional 
phase shift where 5 can change from 0 for a quasi-2D 
cylindrical Fermi surface to 1/8 for a corrugated 3D 
Fermi surface [25]. From the linear extrapolation of 
index plot, the intercept 0.33 (0.20) obtained in Fig. 
4(a) from the SdH oscillation with B along [100] and 
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FIG. 4. (a) Landau index (n) plotted against 1/B (Closed 
circle denote the peak position) from SdH oscillation (From 
the linear extrapolation of fitted line, intercept for B || [100] is 
0.32 and for B jj [021] is 0.20), (b)Landau index n +1/4 plot¬ 
ted against 1/B (Closed circle denote the peak position) from 
dHvA oscillation (Intercept at 0.23 and 0.35 for two parallel 
line) (c) Quantum oscillations of susceptibility Ax(l/B) and 
of resistance AR(1/B) at 10 K temperature, (d)Temperature 
dependence of the relative amplitude of quantum oscillation 
for the 9th Landau level where Ra and Rb are the oscillation 
amplitude in SdH oscillation for B]]]100] and B]][02l] respec¬ 
tively and m corresponds to the oscillation in susceptibility 
.The solid line is a fit to the Lifshitz-Kosevich formula(Equ. 

I). 


([021]). The intercept should be within (+,-)l/8 from 
zero depending on the curvature along z-direction for 
the 3D Dirac fermion system . But we have found an 
additional phase shift S ^ 0.32 (0.2), deviating from the 
expected value. From the index plot of dHvA oscillation, 
we have found two intercepts at 0.23 and 0.35 due to the 
presence of two frequencies. Again in our experiment we 
have used B < 9 T and the Landau index plots are very 
straight, so one cannot expect the phase shifting due to 
the orbital and Zeeman splitting of the spin-degenerate 
conduction band So, we assume that this extra 

phase shift may be due to the contribution from a 
small parabolic curvature and anisotropic ellipsoidal 
nature of Fermi surface away from the Dirac point. This 
result clearly demonstrates a nontrivial Berry’s phase 
associated with 3D Dirac semimetal phase in Cd 3 As 2 . 

Both SdH and dHvA oscillations are appropriate 
tools to study the Fermi surface properties and to 
identify the type of charge carrier. Figure 4(c) shows 
the measured magnetic susceptibility and the oscillating 
part of magnetoresistance as a function of 1/B. The 
oscillation in AR^x is approximately 180° out of phase 


with respect to Ay. As the oscillatory component of 
the conductivity (Aaxx) is out of phase with ARxx, 
i.e., sign(A(Ta:a:)=-sign(ABa;a,), both Acr^x and Ay are 
in phase with each other. Now considering the relation 
Aaxx=fJ'\'m\B^Ax [l^ and taking into account the 
fact that, in- and out-of-phase behavior correspond to 
electrons and holes with fi= +1 and /i = -1 respectively 
, we conclude that the majority charge carrier in this 
system is electron. This is consistent with negative sign 
of thermoelectric power (TEP) of the present samples 
in the temperature range 10-350 K (see supplementary 
material 1^) and the negative sign of reported Hall 
coefficient [^. It may be noted that TEP exhibits a 
linear temperature dependence similar to that observed 
in graphene, which is in quantitative agreement with the 
semiclassical Mott relation [s^. 

The thermal damping of the oscillation amplitude A 
can be described by the temperature dependent part of 
the Lifshitz-Kosevich formula: 

2'K'^kBT/TiWc 

sinh(27r^fcBT /fiUJc) ’ 

where a is a temperature independent constant and 
Wc is the cyclotron frequency. The energy gap 
can be obtained by fitting the relative amplitude of 
SdH and dHvA oscillations with Eq. (1) [Fig. 4(d)]. 
Though for a Dirac system, its linear spectrum implies 
zero rest mass and gives the Landau energy level 
En=VF\/{‘2nheB)^ their cyclotron mass is not zero. 
So the quantum oscillation amplitude can indeed be 
described by the above expression with Wc = eB/m* 
where m/ is the effective cyclotron mass of charge 
carrier and the Fermi velocity up=?tfcp/m* |3ll - l33l| . 
We deduce 1.03x10^(0.96x10®) m/s and m* ^ 
0.045me(0.05me) by analyzing the SdH oscillations with 
field along [021]([100]) direction. Similarly, we have 
fitted the amplitude of 9th Landau Level of dHvA 
oscillations with Eq. (1) and found Mp ~ 1.04x 10®m/s 
and m/~0.044ine for F=53 T. Thus, the calculated 
values of vp and m/ from transport and magnetization 
data are close to each other. These values are also 
close to those reported by different groups [l9l - [2^ .Now. 
using the value of Fermi energy and Fermi velocity, we 
have deduced the Fermi energy ^ 270 meV. So we are 
above the Lifshitz transition point reported in [mil 
By htting the MR oscillation amplitude with magnetic 
field to exp{—2TT'^kBTp)/fiuJc), we deduce the Dingle 
temperature Td= 16.8 K at 2 K. 

In conclusion, we have studied dHvA and SdH oscilla¬ 
tions as two parallel methods to probe the Fermi surface 
of Cd 3 As 2 and the values of Fermi wave vector, Fermi 
velocity, and effective cyclotron mass of charge carrier 
have been calculated from both the techniques. In 
our study, dHvA effect reveals two distinct frequencies 
of oscillations 46 and 53 T along [021], unlike single 
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frequency along [100] and [012], which confirms the ex¬ 
istence of two Fermi surface cross-sections along certain 
directions. Both transport (SdH) and magnetic (dHvA) 
measurements support a non-trivial tt Berry’s phase 
with a phase shift which is the signature of existence of 
a 3D Dirac semimetal phase in Cd 3 As 2 . 

Note added in proof. After submit ting our manuscript, 
we have come to know the results in |34j ] where they have 


found double frequency from SDH oscillation measure¬ 
ment. They have argued that Fermi surface of Cd 3 As 2 
consists two nested ellipsoids beyond the Lifshitz saddle 
point. 
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I. SUPPLEMENTARY MATERIAL FOR 
“PROBING THE FERMI SURFACE OF 3D 
DIRAC SEMIMETAL CD 3 AS 2 THROUGH DE 
HAAS-VAN ALPHEN TECHNIQUE” 

H. CRYSTAL STRUCTURE ANALYSIS 

Phase purity and the structural analysis of the sam¬ 
ples were done by high resolution powder x-ray diffrac¬ 
tion (XRD) technique (Rigaku, TTRAX II) using Cu-Kq, 
radiation. Fig. SI shows the x-ray diffraction pattern 
of powdered sample of Cd 3 As 2 single crystals at room 
temperature. Within the resolution of XRD, we have 
not seen any peak due to impurity phase. Using the Ri- 
etveld profile refinement program of diffraction patterns, 
we have calculated the lattice parameters 0=6=12.64403 
A and c=25.44722 A with space group symmetry Ri/acd. 
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FIG. 5. X-Ray diffraction pattern of powdered single crystals 
of Cd 3 As 2 . In black experimental data, in red the calcu¬ 
lated pattern, in bine the difference and in green the Bragg 
positions. 



FIG. 6. Magnetoresistance of Cd 3 As 2 single crystal with cur¬ 
rent along [012] direction and the magnetic field is applied 
along [100] direction. 


III. MAGNETORESISTANCE MEASUREMENT 

Figure 2 shows the magnetic field dependence of mag¬ 
netoresistance at various temperatures for another piece 
of Cd 3 As 2 single crystal. Similar to Fig. 1(b) in the 
text, the magnetic field was applied along [100] direction 
while the current was along [012]. Though, the value of 
magnetoresistance is small for this sample, the amplitude 
of SdH oscillation is much larger and it is clearly visible 
even at 77 K. 


IV. THERMAL TRANSPORT MEASUREMENT 

Both the resistivity and thermoelectric power (TEP) 
measurements of Cd 3 As 2 single crystals were done by 


standard four-probe technique. Electrical contacts were 
made using conductive silver paste. Measurements were 
done in the temperature range 2-350 K in physical prop¬ 
erty measurement system (Quantum Design). The tem¬ 
perature dependence of TEP is shown in Fig. 3. TEP 



T(K) 


FIG. 7. Temperature dependence of thermoelectric power of 
Gd 3 As 2 single crystal. 

is negative and shows linear T dependence over the mea¬ 
sured temperature range. At room temperature, the 
value of TEP is about 61 /iV/K which is significantly 
large. 

V. MAGNETIZATION MEASUREMENT 



1/B(T') 


FIG. 8. Oscillating part of magnetization Am vs. 1/B plot( 
de Haas van Alphen(dHvA) effect)for magnetic held along 
[100] direction and inset shows oscillation frequency at 55.20 
T after Fast Fourier transform found from dHvA effect. 
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FIG. 9. Oscillating part of magnetization Am vs. XjB plot( 
de Haas van Alphen(dHvA) effect)for magnetic field along 
[012] direction and inset shows the single frequency oscillation 
at 46.5 T after Fast Fourier transform found from dHvA effect. 
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